Measurement of the Neutron Lifetime Using a Proton Trap 
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We report a new measurement of the neutron decay lifetime by the absolute counting of in-beam 
neutrons and their decay protons. Protons were confined in a quasi-Penning trap and counted with 
a silicon detector. The neutron beam fluence was measured by capture in a thin 6 LiF foil detector 
with known absolute efficiency. The combination of these simultaneous measurements gives the 
neutron lifetime: r n = (886.8 ± 1.2[stat] ± 3.2[sys]) s. The systematic uncertainty is dominated by 
uncertainties in the mass of the 6 LiF deposit and the 6 Li(n, t) cross section. This is the most precise 
measurement of the neutron lifetime to date using an in-beam method. 

PACS numbers: 13.30.Ce, 21.10.Tg, 23.40.-s, 26.35.+C 



Precision measurements of neutron beta decay address 
basic questions in particle physics, astrophysics, and cos- 
mology. As the simplest semi-leptonic decay system, the 
free neutron plays a crucial role in understanding the 
physics of the weak interaction and testing the validity 
of the Standard Model p). The current experimental 
uncertainty in the neutron lifetime dominates the uncer- 
tainty in calculating the primordial helium abundance, as 
a function of baryon density, using Big-Bang nucleosyn- 
thesis ■ Figure 1 shows a summary of recent neutron 
lifetime measurements. The four most precise 
confined ultracold neutrons in a bottle; the decay life- 
time is determined by counting the neutrons that remain 
after some elapsed time, with a correction for competing 
neutron loss mechanisms. Two others 0, Q measured 
the absolute specific activity of a beam of cold neutrons 
by counting decay protons. Given the very different sys- 
tematic problems that the two classes of experiments en- 
counter, a more precise measurement of the lifetime using 
the in-beam technique not only reduces the overall un- 
certainty of t„ but also provides a strong check on the 
robustness of the central value. Note that the in-beam 
method measures only the neutron decay branch that 
produces a proton in the final state, while the UCN bot- 
tle method measures all decay modes. With future im- 
provements in precision it may become possible to place 
important limits on exotic, nonbaryonic branches of neu- 
tron decay by comparing the lifetimes measured by the 
in-beam and UCN bottle methods. 

Figure 2 shows a sketch of the experimental method, 
ft is based on the method first proposed by Byrne, et 
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FIG. I: Figure I: A comparison of recent neutron lifetime 
measurements. The shaded band is ±1 standard deviation of 
the weighted average, including our result. 



al, and is described in greater detail in previous publica- 
tions ^3] • A proton trap of length L intercepts the entire 
width of the neutron beam. Neutron decay is observed 
by trapping and counting decay protons within the trap 
with an efficiency e p . The neutron beam is characterized 
by a velocity dependent fluence rate I(v). The mean 
number of neutrons inside the trap volume is 



N n = L [ dal(v)-, 

J A V 



(1) 



where A is the beam cross sectional area, and so the rate 
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FIG. 2: Figure 2: Schematic representation of the experimen- 
tal method (not to scale). 



N p at which decay protons are detected is 



.V„ =^ f dal(v)-. 



(2) 



After leaving the trap, the neutron beam passes through 
a thin foil of 6 LiF. The probability for absorbing a neu- 
tron in the foil through the 6 Li(n,i) 4 He reaction is in- 
versely proportional to the neutron velocity v. The 
reaction products, alphas or tritons, are counted by a 
set of four silicon surface barrier detectors in a well- 
characterized geometry. We define the efficiency for the 
neutron detector, e , as the ratio of the reaction prod- 
uct rate to the neutron rate incident on the deposit for 
neutrons with thermal velocity v — 2200 m/s. The cor- 
responding efficiency for neutrons of other velocities is 
£ Vo/v. Therefore, the net reaction product count rate 
N n is 



N a = e v I dal(v)-. 

' A V 



(3) 



The integrals in Eq. J5J and Eq. are identical; the 
velocity dependence of the neutron detector efficiency 
compensates for the fact that the faster neutrons in the 
beam spend less time in the decay volume. This cancel- 
lation is exact given two assumptions: (1) the neutron 
absorption efficiency in the 6 LiF target is exactly pro- 
portional to 1 jv and (2) the neutron beam intensity and 
its velocity dependence do not change between the trap 
and the target. The deviation from the l/v law in the 
6 Li(n, t) 4 He cross section at thermal and subthermal en- 
ergies has been shown to be less than 0.01% [llj, and 
changes in the neutron beam due to decay in flight and 
residual gas interaction are less than 0.001%. Our only 
non-negligible correction to Eqs. (J2J and (fJJ) is for the 
finite thickness of the 6 LiF foil (+5.4 s). Thus we obtain 
the neutron lifetime r„ from the experimental quantities 
N a /Np, e D , e p , and L. 

The experiment was performed at the cold neutron 
beamline NG6 at the National Institute of Standards and 
Technology (NIST) Center for Neutron Research. Cold 
neutrons were transported 68 m from the cold source to 
the experiment by a straight 58 Ni-coated neutron guide. 



Immediately after exiting a local guide shutter, the neu- 
tron beam passed through 20 cm to 25 cm of cooled single 
crystal bismuth to attenuate fast neutrons and gamma 
rays that contribute to the background signal. The beam 
was collimated by two 6 LiF apertures. The diameter of 
the first aperture (CI) was varied from 32 mm to 51 mm 
as a systematic check on the effect the beam diameter on 
the measured lifetime. The second (C2) had a diameter 
of 8.4 mm and was not changed. In between these two 
beam-defining apertures were several 6 LiF beam scrap- 
ers, which removed scattered and very divergent neu- 
trons. After C2, the beam transited aim section of 
pre-guide and then entered the vacuum system through 
a 7.5 mm inner diameter quartz guide tube. After passing 
through the trap, it traveled 83 cm to the 6 LiF deposit. 
The beam exited the vacuum system through a silicon 
window and stopped in a thick 6 LiF beam dump. 

The proton trap was a quasi-Penning trap, composed 
of sixteen annular electrodes, each 18.6 mm long with 
an inner diameter of 26.0 mm, cut from fused quartz 
and coated with a thin layer of gold. Adjacent segments 
were separated by 3 mm-thick insulating spacers of un- 
coated fused quartz. The dimensions of each electrode 
and spacer were measured to a precision of ±5 /im using 
a coordinate measuring machine at NIST. Changes in the 
dimension due to thermal contraction are below the 10 
level for fused quartz. The trap resided in a 4.6 T mag- 
netic field and the vacuum in the trap was maintained 
below 10~ 9 mbar. 

In trapping mode, the three upstream electrodes (the 
"door") were held at +800 V, and a variable number 
of adjacent electrodes (the "trap") were held at ground 
potential. The subsequent three adjacent electrodes (the 
"mirror" ) were held at +800 V. We varied the trap length 
from 3 to 10 grounded electrodes. When a neutron de- 
cayed inside the trap, the decay proton was trapped radi- 
ally by the magnetic field and axially by the electrostatic 
potential in the door and mirror. After some trapping pe- 
riod, typically 10 ms, the trapped protons were counted. 
The door was "opened", i.e. the door electrodes were 
lowered to ground potential, and a small ramped poten- 
tial was applied to the trap electrodes to assist the slower 
protons out the door. The protons were then guided by 
a 9.5° bend in the magnetic field to the proton detector 
held at a high negative potential (-27.5 kV to -32.5 kV). 
After the door was open for 76 /is, a time sufficient to 
allow all protons to exit the trap, the mirror was lowered 
to ground potential. This prevented negatively charged 
particles, which may contribute to instability, from ac- 
cumulating in any portion of the trap. That state was 
maintained for 33 /is, after which the door and mirror 
electrodes were raised again to +800 V and another trap- 
ping cycle began. Since the detector needed to be enabled 
only during extraction, the counting background was re- 
duced by the ratio of the trapping time to the extraction 
time (typically a factor of 125). Figure 3 shows a plot of 
proton detection time for a typical run. 

Protons that were born in the trap (grounded electrode 
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FIG. 3: Figure 3: A typical plot of proton detection time after 
gating on the detector (t — 0). Regions I and III are used to 
subtract background from the peak region II. 



region) were trapped with 100% efficiency. However pro- 
tons that were born near the door and mirror (the "end 
regions"), where the electrostatic potential is elevated, 
were not all trapped. A proton born in the end region 
was trapped if its initial (at birth) sum of electrostatic 
potential energy and axial kinetic energy was less than 
the maximum end potential. This complication caused 
the effective length L of the trap to be difficult to de- 
termine precisely. It is for this reason that we varied 
the trap length. The shape of the electrostatic potential 
near the door and mirror was the same for all traps with 
3-10 grounded electrodes, so the effective length of the 
end regions, while unknown, was in principle constant. 
The length of the trap can then be written L = nl + L cn d 
where n is the number of grounded electrodes and I is the 
physical length of one electrode plus an adjacent spacer. 
Lend is an effective length of the two end regions; it is pro- 
portional to the physical length of the end regions and 
the probability that protons born there will be trapped. 
From Eq. and © we see that the ratio of proton 
counting rate to alpha counting rate is then 



N, 



p =T -i 



(nl + Lend)- 



(4) 



We fit N p /N a as a function of n to a straight line and 
determine r n from the slope, so there is no need to know 
the value of Lend, provided that it was the same for all 
trap lengths. 

Because of the symmetry in the Penning trap's de- 
sign, Lend was approximately equal for all trap lengths. 
There were two trap-length-dependent effects that broke 
the symmetry: the gradient in the axial magnetic field, 
and the divergence of the neutron beam. Each of these 
effects caused L en d to vary slightly with trap length. A 
detailed Monte Carlo simulation of the experiment, based 
on the measured and calculated magnetic and electric 
field inside the trap, was developed in order to correct 
for these trap nonlinearities. It gave a trap-length de- 
pendent correction that lowered the lifetime by 5.3 s. 



Surface barrier (SB) and passivated ion-implanted pla- 
nar silicon (PIPS) detectors were used at various times 
for counting the protons. The proton detectors were large 
enough so that all protons produced by neutron decay in 
the collimated beam, defined by CI and C2, would strike 
the 19.7-mm diameter active region after the trap was 
opened. The detector was optically aligned to the mag- 
netic field axis, and the alignment was verified by scan- 
ning with a low energy electron source at the trap's center 
and with actual neutron decay protons. When a proton 
hit the active region of the detector, the efficiency for pro- 
ton detection was less than unity because: 1) a proton 
can lose so much energy in the inactive surface dead layer 
that it deposits insufficient energy in the active region to 
be detected; 2) a proton can be stopped within the dead 
layer and never reach the active region; 3) a proton can 
backscatter from the dead layer or the active region and 
deposit insufficient energy. In the first two cases, the 
proton will definitely not be detected. We denote the 
fraction of protons lost in this manner by /Lost ■ In the 
last case, there is some probability that the backscattered 
proton will be reflected back to the detector by the elec- 
tric field and subsequently be detected. We denote the 
fraction of protons that backscatter but still have some 
chance for detection by /bsc- 

To determine the proton detection efficiency, we ran 
the experiment with a variety of detectors with different 
dead layer thicknesses and different acceleration poten- 
tials. The fraction of protons lost /Lost and the fraction 
that backscatter /b sc were calculated using the SRIM 
2000 Monte Carlo program 0. We found that / Lost 
varied from 4.0(3) x 10~ 5 to 8.0(6) x 10~ 3 , and / Bsc from 
1.83(13) x 10~ 3 to 2.37(17) x 10~ 2 , depending on the de- 
tector and acceleration potential. The proton counting 
rate N p for each run was multiplied by 1 + /Lost to correct 
for lost protons. The correction for backscattered protons 
was not as simple because of the unknown probability for 
a backscattered proton to return and be detected, so we 
made an extrapolation of the measured neutron lifetime 
to zero backscatter fraction (see Figure 4). 

The neutron detector target was a thin (0.34 mm), 
50-mm-diameter single crystal wafer of silicon coated 
with a 38 mm diameter deposit of 6 LiF, fabricated at 
the Institute for Reference Materials and Measurements 
(IRMM) in Geel, Belgium. The manufacture of de- 
posits and characterization of the 6 LiF areal density 
were exhaustively detailed in measurements performed 
over several years 01 ■ The average areal density was 
p = (39.30 ± 0.10) /ig/cm 2 . The a particles and tritons 
produced by the the 6 Li(n, i) 4 He reaction were detected 
by four surface barrier detectors, each with a well-defined 
and carefully measured solid angle. The geometry was 
chosen to make the solid angle subtended by the detec- 
tors insensitive to first order in the source position. The 
parameter en, gives the ratio of detected alphas/tritons to 
incident thermal neutrons. It was calculated using 



= irr I I &(x,y)p(x,y)6(x,y)dxdy, ( 5 ) 
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FIG. 4: Figure 4: A linear fit of the measured neutron lifetime 
at varying values of the detector backscattering fraction. The 
extrapolation to zero backscattering gives the free neutron 
lifetime. 

where <7o is the cross section at thermal (vq = 2200 m/s) 
velocity, Q(x,y) is the detector solid angle, p(x,y) is the 
areal mass density of the deposit, and 0(x, y) is the areal 
distribution of the neutron intensity on the target. The 
6 Li thermal cross section is (941.0 ± 1.3) b 01- It i s lm ~ 
portant to note that we take the Evaluated Nuclear Data 
Files (ENDF/B-6) la uncertainty from the evaluation, 
not the expanded uncertainty, to be the most appropri- 
ate for use with this precision experiment. The neutron 
detector solid angle was measured in two independent 
ways: mechanical contact metrology and calibration with 
239 Pu alpha source of known absolute activity. These two 
methods agreed to within 0.1 %. 

The neutron beam intensity distribution was measured 
at both ends of the trap, and at the neutron monitor, 
using the Dy foil imaging technique |15|. This method 
provides almost 5 decades of dynamic range, making it 
ideal for sensitive neutron measurements. These quan- 
titative images were used to determine the beam distri- 
bution function 6(x,y) in Eq. [SJ They were also used 
to measure the beam "halo" , the scattered neutrons that 
lie outside the collimated neutron beam. Trapped decay 
protons that originate from the halo may have trajecto- 
ries that lie outside the active radius of the proton detec- 
tor so they are not counted, while these halo neutrons are 
counted by the neutron counter. We found that 0.1 % of 
the neutron intensity is in the halo outside the effective 
radius of the proton detector. 

Proton and neutron counting data were collected for 
13 run series, each with a different proton detector and 
acceleration potential. The corrected value of the neu- 
tron lifetime for each series was calculated and plotted 
vs. backscattering fraction, as shown in Figure 4. A lin- 
ear extrapolation to zero backscattering gave a result of 
t„ = (886.8 ± 1.2[stat] ± 3.2[sys]) s. A summary of all 
corrections and uncertainties is given in Table [I] This 
result would be improved by an independent absolute 
calibration of the neutron counter, which would signifi- 
cantly reduce the two largest systematic uncertainties, in 
the 6 LiF foil density and 6 Li cross section. A cryogenic 
neutron radiometer that promises to be capable of such 



TABLE I: Summary of systematic corrections and uncertain- 
ties (in seconds) for the measured neutron lifetime. 



Source 


Correction 


Uncertainty 


6 LiF foil areal density 




2.2 


6 Li cross section 




1.2 


Neutron detector solid angle 


+1.5 


1.0 


Neutron beam halo 


-1.0 


1.0 


LiF target thickness 


+5.4 


0.8 


Trap nonlinearity 


-5.3 


0.8 


Neutron losses in Si wafer 


+1.3 


0.5 


6 LiF distribution in deposit 


-1.7 


0.1 


Proton backscatter calc. 




0.4 


Proton counting statistics 




1.2 


Neutron counting statistics 




0.1 


Total 


+0.2 


3.4 



a calibration at the 0.1% level has recently been demon- 
strated 0|, and we are pursuing this method further. 
With this calibration, we expect that our experiment will 
ultimately achieve an uncertainty of less than 2 seconds 
in the neutron lifetime. 
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